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Abstract. r/ Carinae is considered to be a massive colliding wind binary system with 
a highly eccentric ( e ~ 0.9), 5.54-yr orbit. However, the companion star continues 
to evade direct detection as the primary dwarfs its emission at most wavelengths. Us- 
ing three-dimensional (3-D) SPH simulations of r\ Car’s colliding winds and radiative 
transfer codes, we are able to compute synthetic observables across multiple wavebands 
for comparison to the observations. The models show that the presence of a companion 
star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line 
profiles, as well as the ground-based photometric monitoring. Here, we focus on the 
Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in 
the dense primary wind, allowing increased escape of radiation from the hotter/deeper 
layers of the primary’s extended wind photosphere. The results have important impli- 
cations for interpretations of tj Car’s observables at multiple wavelengths. 


1. Introduction 

The 5.54-year periodicity observed in numerous spectral lines, as well as the 1R, Vi- 

sual, and X-ray fluxes, strongly suggests r\ Carinae is a binary system (Whitelock et al. 

2004; Damineli et al. 2008a, b; Femandez-Lajus et al. 2010; Corcoran 201 1). Recently, 
Okazaki et al. (2008, hereafter 008) modeled t] Car’s RATE X-ray light curve using 
a three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulation of the 

binary' wind-wind collision. A key point of 008 is that the fast wind of the secondary- 
star. ij B , canes a cavity in the dense wind of the primary, q A , allowing X-rays that 
would othenvise be absorbed to escape into our line-of-sight. This immediately sug- 
gests that if the primary wind is sufficiently optically thick in the UV, Optical, or IR 
waveband, the low-density secondary wind may likewise carve or “bore” a cavity or 
“hole” in the associated wind photosphere, allowing increased escape of radiation from 
the hotter/deeper layers. 
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2. The Bore-Hole Effect 


Any bore-hole effect should depend on (1) how close the wind cavity carved by ijb 
gets to t] A and (2) the apparent size of // j’s photosphere. For a bore-hole effect to 
occur, the optically thick primary wind photosphere must extend at least as far as the 
separation distance between the primary and the apex of the wind-wind interaction 
region (WWIR), i.e. if at some point (2) > (1), there is a bore hole. 

The separation between tja and the apex of the WWIR is determined by the dy- 
namics of the colliding winds and the orbit. The absolute minimum distance that rj A 's 
photosphere must extend is the separation between i] A and WWIR apex at periastron: 


a{ 1 — e) 

= 1 + 1/V^’ 


(1) 


where e is the eccentricity (0.9), a is the orbital semi-major axis length (15.4 AU), and 
rj is the momentum ratio of rj/ s wind to //g’s wind. For the parameters of 008, tj ~ 4.2 
and R m j n ~ 1 AU. Using the primary mass-loss rate M = 10 " 3 M Q yr" 1 suggested by 
Hillier et al. (2001, 2006, hereafter H01 and H06), n « 16.67 and R mia ~ 1.24 AU. 

The apparent size of rj/s wind photosphere is determined by radiative transfer 
effects. While there are many ways to define this apparent size, a simple way is to use 
the radial photospheric radius 7? p h 0 t at which the radial optical depth r = 1 : 


_ kM 

Rphot ~~ 4/rv ’ 


( 2 ) 


where k is the opacity (units of cm 2 g" 1 , assumed constant for simplicity), M is the 
mass-loss rate of q A , and v is the terminal speed of rj A ’s wind. Assuming k = 1 cm 2 g" 1 
and mass-loss rates of M = 2.5x 10 “ 4 M Q yr" 1 and M = 10" 3 M G yr" 1 , R p h 0 t ~ 1-67 AU 
and 6.71 AU, respectively. Thus, even for low values of k of order unity, a bore-hole 
effect should occur in the rj Car system (at least at periastron). 

In reality, any bore-hole effect in r] Car will be wavelength dependent and occur 
over a relatively short period of time around periastron. The wavelength dependence 
is due to the opacity, which determines the size of rj/s photosphere. Figure 1 illus- 
trates this, showing an overlaid plot of the normalized intensity as a function of impact 
parameter from the center of tja in each of the wavebands B, V, R, 1, J, H, K, L. 

Since the separation between tj a and the apex of the WWIR varies with orbital 
phase (p, the bore-hole effect is also orbital-phase dependent. At apastron, the separa- 
tion is ~ 23.6 AU, much larger than jj a in either the K- or B- band. Because of the large 
orbital eccentricity, this separation remains large for most of the orbit, > 10 AU for (f> 
between 0.5 and 0.925. However, between <p w 0.925 and 1.0 (periastron), the separa- 
tion drops rapidly, from ~ 10 AU to 1.24 AU. At tp ~ 0.97, the separation is ~ 5 AU, 
which is about the same size as t] A in the R'-band. Therefore, for the wavebands con- 
sidered here, a bore-hole effect should occur in ?/ Car at phases 0.97 < ,<p < 1 .03. 

Because of the different photospheric radii in each waveband, the wind cavity will 
bore different relative distances into r] A at a given d> in a particular waveband. Thus, at 
cb = 0.98, one would expect to see a significant bore-hole effect in the Al-band, but only 
a relatively minor effect, if any, in the S-band. Moreover, because of the differences 
in relative flux probed between wavebands, one might also expect different relative 
changes in brightness before, during, and after periastron. 
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Figure 1. Curves showing normalized intensity vs. impact parameter in each of 
the wavebands (from innermost to outermost) B, V, R, /, J, H, K, L, computed using 
the HOI, H06 CMFGEN model of t]a- 



Figure 2. 3-D renderings of the bore-hole effect at orbital phases of apastron (left), 

5 days before periastron (middle), and 155 days after periastron (right). The x and 
y axes are the major and minor axes, respectively, and z is the orbital axis X to 
the orbital plane. Lengths are in semi-major axes and grayscale indicates surface 
brightness. See Madura (2010) for modeling details. 


To investigate the bore-hole effect, we use 3-D SPH simulations of ?/ Car’s col- 
liding winds, together with a modified version of the visualization program SPLASH 
(Price 2007), to generate renderings of the surface brightness of r\ A in each of the var- 
ious wavebands B through L as a function of orbital phase. Opacities and primary star 
source functions as a function of radius in each waveband were taken from the spheri- 
cally symmetric model of r] A by HOI, H06. Details can be found in Madura (2010). 

Our models reveal that when tjb is at or near apastron, the WWIR is indeed too far 
from rj A and there is no bore hole (Fig. 2, left panel). As j]b moves closer to rj A during 
its orbit, the WWIR gradually penetrates into the primary photosphere, creating a bore- 
hole effect that increases up until periastron passage (Fig. 2, middle panel), at which 
point i]b quickly wraps around the back side of t]a, the bore hole briefly vanishing as it 
faces away from the observer. After periastron, the bore hole reappears on the opposite 
side of tj A (Fig. 2, right panel) and then slowly fades as t]b moves back toward apastron. 

Using the above 3-D model, we also generate synthetic photometric light curves in 
each waveband for comparison to the available ground-based photometry (Whitelock et al. 
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Figure 3. Synthetic light curves from the bore-hole model (circles) com- 
pared to the photometric observations of q Car’s 2003.5 event (squares, 
Femandez-Lajus et al. 2010) in the B (left) and K (right) bands . 


2004; Femandez-Lajus et al. 2010). We find that our bore-hole model reproduces the 
steep rise and drop observed in each waveband before periastron, and gives roughly 
the same peak-to-peak change in magnitude and duration of the ‘eclipse-like’ events 
seen every 5.54-years (see Figure 3 and Madura 2010). We thus find that t)B and the 
WWIR have a profound influence on q/s extended photosphere and the photometric 
observables. The bore-hole effect likely explains q Car’s periodic eclipse-like events. 


3. Effects of a Bore Hole on // Car’s Spectrum and Interferometric Observables 

Using the 2-D code of Busche & Hillier (2005), we have also investigated how a bore- 
hole effect could alter the observed optical and UV spectra of // Car. Details on the 
modeling approach can be found in Groh (201 1 ) and references therein. 

Assuming a single-star scenario, HO 1 obtained a reasonable fit to the observed HST 
optical spectrum obtained just after periastron, <p ~ 0.05. The model spectrum from 
H01 reproduces well the emission line profiles of H, Fe II, and N II lines. However, 
compared to the observations, the H01 model overestimates the amount of P-Cygni 
absorption. The comparison is even worse as one moves toward apastron, when the 
observations show little or no P-Cygni absorption in H and Fe II lines. 

Using our 2-D radiative transfer model, which takes into account the cavity in 
the wind of q A caused by tjb, we computed the synthetic optical spectrum of q Car at 
4> ~ 0.6 (see figure 1 of Groh 2011). We assume the same parameters for q A as H01 and 
a standard geometry of the WWIR at apastron as predicted by the 3-D SPH simulation. 
For a viewing angle with inclination i = 41° and longitude of periastron co = 270°, the 
2-D model produces a much better fit to the P-Cygni absorption line profiles of H and 
Fe II lines than the l-D HOI model, while still fitting the emission line profiles. The 
improved fit to the P-Cygni absorption is due to the cavity in the wind of q A , which 
reduces the H and Fe II optical depths in line-of-sight to the primary. We find similar 
results when modeling q Car's observed UV spectrum. 

A preliminary study of synthetic Ha line profiles generated using our 2-D cav- 
ity model of q A also reveals broader, deeper P-Cygni absoiption near the stellar poles, 
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Figure 4, Synthetic line profiles of Her generated using the 2-D model discussed 
in the text, emphasizing the P-Cygni absorption at three stellar latitudes for ij A . 


compared to the equator (Figure 4). A bore-hole or wind-cavity model may therefore 
provide an alternative explanation for HST long-slit spectral observations of Ha line 
profiles as a function of stellar latitude that have been historically interpreted as evi- 
dence that the wind of tj A is latitude dependent (Smith et al. 2003). 

Finally, using a similar 2-D model, Groh et al. (2010) find that the density struc- 
ture of j] A ’s wind can be sufficiently disturbed by i] R and the associated wind cavity, thus 
mimicking the effects of fast rotation in the interferometric observables. Groh et al. 
(20 1 0) further show that even if q A is a rapid rotator, models of the interferometric data 
are not unique, and both prolate- and oblate-wind models can reproduce current inter- 
ferometric observations. These prolate- and oblate-wind models additionally suggest 
that the rotation axis of rj A would not be aligned with the Homunculus polar axis. 


4. Conclusions 

Our conclusions are rather straightforward; based on the results of our 3-D hydrody- 
namical simulations and multi-D radiative transfer modelling, we find that the sec- 
ondary star i]b and WW1R significantly affect the extended wind photosphere of r] A , 
and thus multi-wavelength observations of the // Car system. 

Acknowledgments. We thank Nathan Smith for discussions that led to this work. 
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